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ABSTRACT 

All sufficiently massive clusters of galaxies are expected to be surrounded by strong accretion shocks, 
where protons can be accelerated to ~ 10 18 -10 19 eV under plausible conditions. Such protons inter- 
act with the cosmic microwave background and efficiently produce very high energy electron-positron 
pairs, which then radiate synchrotron and inverse Compton emission, peaking respectively at hard 
X-ray and TeV gamma-ray energies. Characterized by hard spectra (photon indices ~ 1.5) and spa- 
tial distribution tracing the accretion shock, these can dominate over other nonthermal components 
depending on the shock magnetic field. HESS and other Cerenkov telescopes may detect the TeV emis- 
sion from nearby clusters, notwithstanding its extended nature. The hard X-rays may be observable 
by future imaging facilities such as NeXT, and possibly also by ASTRO-E2/HXD. Such detections 
will not only provide a clear signature of ultra-high energy proton acceleration, but also an important 
probe of the accretion shock itself as well as magnetic fields in the outermost regions of clusters. 
Subject headings: acceleration of particles — cosmic rays — radiation mechanisms: non-thermal — 
gamma rays: theory — galaxies: clusters: general — X-rays: galaxies: clusters 



1. INTRODUCTION 

The standard cosmological picture of hierarchical 
structure formation dictates that virialized objects such 
as clusters of galaxies continue to grow by accreting dark 
matter and baryonic gas. It is generically expected that 
strong accretion shocks occur around the virial radii of 
all sufficiently massive clusters, where inf ailing gas of rel- 
atively low temperature is heated up to the virial value 
(e.g. Voit 2004). Here "accretion" signifies not only in- 
fall of diffuse gas but also minor mergers, i.e. merging of 
virialized objects with relatively small masses compared 
to the cluster. Such shocks are potentially interesting 
sites of high energy particle acceleration and nonther- 
mal emission, possessing high Mach numbers and capa- 
ble of generating particles with hard spectra at high ef- 
ficiency (Ryu et al. 2003; Gabici & Blasi 2003; Inoue & 
Nagashima 2005). This is in contrast to major merger 
shocks, which are likely to play a lesser role for high en- 
ergy phenomena on account of their low Mach numbers 
(c.f. Totani & Kitayama 2000; Takizawa & Naito 2000; 
Berrington & Dermer 2003). 

Primary electrons, i.e. those accelerated directly by 
the shocks, give rise to radio synchrotron emission as 
well as inverse Compton (IC) emission by upscattering 
the cosmic microwave background (CMB) to UV, hard 
X-ray and gamma-ray energies (e.g. Hwang 1997; Enfilin 
& Biermann 1998; Sarazin 1999; Atoyan & Volk 2000; 
Waxman & Loeb 2000; Petrosian 2001). Such emis- 
sion should spatially trace the morphology of their in- 
jection sites, which are shell-like regions of Mpc radii 
for accretion shocks (e.g. Miniati 2003). On the other 
hand, protons interact with the intracluster gas and gen- 
erate ir° decay gamma-ray emission (Volk, Aharonian & 
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Breitschwerdt 1996; Berezinsky, Blasi & Ptuskin 1997; 
Pfrommer & Enfilin 2004), as well as synchrotron and IC 
emission by secondary electrons and positrons from 7r 
decay (Blasi & Colafrancesco 1999), all of which should 
follow the centrally-concentrated spatial distribution of 
thermal gas. 

The current observational evidence for such kinds of 
nonthermal emission is still ambiguous. Various types of 
diffuse radio emission in clusters are known (Kempncr 
et al. 2003; Feretti, Burinaga & Enfilin 2004); however, 
except for a few cases, they cannot be easily explained as 
due to primary electrons nor secondary electrons from p- 
p interactions (Brunetti 2004). More controversially, dif- 
fuse excess emission over the known thermal component 
has been reported for a few clusters in the extreme UV 
band (e.g. Bowyer et al. 2004) and in hard X-rays (e.g. 
Rephaeli & Gruber 2002; Fusco-Femiano et al. 2004; 
Henriksen & Hudson 2004), which has been interpreted 
as IC emission from either primary or secondary elec- 
trons. There is no clear case yet for gamma-ray emission 
from any cluster (Reimer et al. 2003). 

Cluster accretion shocks have also been proposed 
as candidate sources of ultra-high energy (UHE) cos- 
mic rays (Norman, Melrose & Achterberg 1995; Kang, 
Rachen & Biermann 1997). However, despite being vi- 
able on energetic grounds, realistic evaluations of the 
maximum proton energy seem to fall short of 10 20 eV 
by 1 to 2 orders of magnitude (Ostrowski & Siemieniec- 
Ozieblo 2000; Jones 2004). 

In this letter, we point out that protons plausibly accel- 
erated to 10 18 -10 19 eV in cluster accretion shocks should 
induce important radiative signatures that have not been 
investigated so far: synchrotron and IC emission from 
secondary electron-positron pairs produced in p-7 inter- 
actions with the CMB. Aharonian (2002) and Rordorf, 
Grasso & Dolag (2004) have discussed these processes by 
assuming that some source, e.g. a powerful radio galaxy, 
resides inside a cluster and injects UHE protons into its 
surroundings. Here such mechanisms are studied in the 
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context of protons originating in accretion shocks that 
may be generic to all massive clusters. We adopt the 
WMAP cosmological parameters (Spergel et al. 2003): 
Cl m = 0.27, Q A = 0.73, h = 0.71 and Q b = 0.044. 

2. MODEL 

For concreteness, we consider a fiducial cluster of mass 
M = 2 x 10 15 Af Q located at distance D = 100 Mpc, 
similar to the Coma cluster. Spherical symmetry and 
steady state conditions are assumed for simplicity. The 
accretion shock radius R s is taken to be the virial radius, 
R s ~ 3.2 Mpc (Voit 2004). Since gas falls in at a velocity 
equal to the circular velocity V c = {GM/R S ) 1/2 ~ 1600 
km/s and comes to rest downstream of the accretion 
shock as the shock expands outwards, the shock veloc- 
ity relative to the upstream gas is 14 — (4/3)14 — 2200 
km/s. The magnetic field at the shock B s is uncertain, as 
Faraday rotation measurements do not give strong con- 
straints at radii > 1 Mpc (Clarke, Kronberg & Bohringer 
2001, Johnston-Hollitt & Ekers 2004). Some observa- 
tions of diffuse synchrotron emission in intergalactic re- 
gions outside but near clusters suggest field strengths 
— 0.1-1/zG (Kronberg 2004). We thus choose B s to be 
a parameter with 0.1-1/zG as a reasonable range, corre- 
sponding to 0.05-5 % of the gas thermal energy at R s for 
our fiducial cluster, when an isothermal density profile 
and gas fraction f g = Qb/Q rn = 0.16 are assumed. 

We first evaluate the maximum accelerated proton en- 
ergy as in Kang et al. (1997). The timescale t acc for 
accelerating a proton up to energy E p in strong shocks 
is (Blandford & Eichler 1987) 

_ 20 r\r g c 

tacc~ 3 

where r\ characterizes the mean free path for scattering 
by magnetic disturbances, normalized by the particle gy- 
roradius r g = E p /eB s . Observations of nonthermal emis- 
sion from supernova remnant (SNR) shocks indicate that 
rj is close to unity, i.e. the Bohm diffusion limit (e.g. 
Volk, Berezhko & Tsenofontov 2005). Although it is not 
guaranteed that this condition is also realized in cluster 
accretion shocks, we will assume rj = 1 in what follows. 

On the timescale of Eq. Q , energy losses by photopair 
and photopion interactions with the CMB become signif- 
icant. We use simple fitting formulae for the timescales 
of each, which agree with more accurate calculations (e.g. 
Yoshida & Teshima 1993) to within ~ 20 % in the per- 
tinent energy ranges: t pl = max[to exp(_E c /£ , p ), t\], with 
t = 3.2 x 10 9 yr, t x = 4.6 x 10 9 yr and £ c = 3x 10 18 eV 
for photopair losses, and to = 2.65 x 10 7 yr, t\ — 4.5 x 10 7 
yr and E c = 2.65 x 10 20 eV for photopion losses (see also 
Waxman 1995). In addition, we consider the dynami- 
cal timescale td yn — Rs/Vc — 1-9 Gyr to be the life- 
time of the accretion shock, limiting the entire process. 
On timescales much longer than this, the conditions may 
change so appreciably as to invalidate our steady state 
approximation (e.g. the cluster may undergo a major 
merger). Note that the timescale for adiabatic losses 
in the expanding accretion shock will always be larger, 
~ Stdym as the shock expansion velocity is ~ 14/3. 



The maximum proton energy i? max is obtained by 
equating t acc and min[£ P7 , td yn ]- With B s = 0.1-1//G, 
£ max ~ 4.5 x 10 18 eV (B s /l/iG) (i dyn /1.9Gyr), always 
limited by the shock lifetime. Although this fails to reach 
the realm of the observed highest energy cosmic rays, it 
is noteworthy that t pl is within an order of magnitude of 
tdyn in the range E p ~ 10 18 -10 19 eV. This implies that 
during the process, a significant fraction of the proton 
energy can be channeled into secondary pairs and conse- 
quent radiation. 

The protons are assumed to form an energy distribu- 
tion with spectral index 2, appropriate for strong shocks 
in the test particle limit, and an exponential cutoff at 
E max , i.e. tx E p 2 exp(-E p /E max ). This is injected 
downstream of the shock for a duration £j n j = td yn , at 
a constant rate normalized by the proton luminosity L p . 
We posit L p to be a fraction f p = 0.1 of L acc , the kinetic 
energy flux through strong accretion shocks, which can 
be expressed 

Lacc — fgGMM/ R s = /,,/:,.. I /l/' //, 

-2.9 x 10 46 erg/s fg facc ( ^- (2) 

8/ 0.16 0.1 \2x 10 15 Af Q y V ' 

where M = f acc Vc /G is the mass accretion rate (White 
1994), and f acc ~ 0.1 is a factor that was estimated by 
Keshet et al. (2004) using numerical simulations. The 
numerical value above is for our chosen cosmological pa- 
rameters, and is consistent with the simulation results of 
Miniati (2002) as well as the analytic estimates of Gabici 
& Blasi (2004a). 

Given the proton spectrum and injection rate, we fol- 
low the method described by Aharonian (2002) to com- 
pute secondary particle production from proton inter- 
actions and the resultant emission. The emission re- 
gion is taken to be a uniform sphere of radius R s for 
simplicity, although this should actually be a shell be- 
hind the shock. The kinetic equations for the energy 
distributions of protons and secondary pairs are solved 
self-consistently with the emitted spectra, including the 
processes of p-7 pair production, p-7 pion production, 
synchrotron and IC radiation from all pairs, and proton 
synchrotron radiation. Gamma-ray and pair production 
from p-p interactions can also be calculated for a given 
ambient gas density; we take this to be n = 10 -6 cm~ 3 , 
appropriate for the low density regions around R s . The 
diffusion coefficient for protons inside the emission region 
is considered to be similar to that at the shock, imply- 
ing that diffusive escape is unimportant up to energies 
E esc ~ 10_ 19 eV 77" 1 (S s /MG)(i? s /3Mpc) 2 (t m ,/2Gyr)- 1 . 
Intergalactic gamma-ray absorption is incorporated us- 
ing the diffuse extragalactic infrared background (IRB) 
model discussed in Aharonian (2001). 

3. RESULTS 

Fig^shows the broadband spectra of emission induced 
by protons in the accretion shock of our fiducial cluster, 
for the cases of B s = 0.1, 0.3 and l^G. These respec- 
tively imply S max = 4.5 x 10 17 , 1.4 x 10 18 and 4.5 x 10 18 
eV, close to or above the CMB photopair threshold but 
below the photopion threshold. The pairs are injected 
at very high energies — (2m e /m p )E p — 10 15 eV and 
then rapidly lose energy via synchrotron and IC cool- 
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ing, developing an energy distribution oc E~ 2 and emit- 
ting hard spectra with photon indices T ~ 1.5 (Aharo- 
nian 2002). Spectral peaks occur at ~ 1-100 keV for 
synchrotron and at ~ 10-100 TeV for IC, while inter- 
galactic absorption by the IRB imposes a further cutoff 
above ~ 10 TeV in the observed flux. Note the large 
influence of B s , which affects both the synchrotron cool- 
ing rate and the peak energy as well as the amount of 
pair injection by changing the fraction of protons above 
threshold through -E max - Lower B s causes lower E msx 
and less pair injection but also reduces the synchrotron 
cooling, leading overall to increased IC and decreased 
synchrotron luminosities; for higher B s , vice- versa. The 
relative luminosities together with the peak frequency lo- 
cations should allow B s to be reliably determined from 
observations. For our range of B Sl the emitted power 
is mostly dominated by IC gamma-rays and can be as 
large as ~ 10 44 erg/s, despite Klein-Nishina suppression 
being effective near the peak. This testifies that protons 
above threshold are radiatively efficient in the sense that 
a major portion of their energy is converted to radiative 
pairs. With n = 10 _6 cm~ 3 , ir° gamma-rays and sec- 
ondary pair radiation from p-p interactions only make 
a miniscule contribution, although they should be more 
pronounced from the central, higher density regions of 
the cluster, which is not considered here. 

Owing to the hard spectra and high radiative effi- 
ciency from protons, the p-7 pair (hereafter PGP) emis- 
sion may dominate the hard X-ray and TeV gamma- 
ray bands over other components from primary elec- 
trons or p-p interactions, which should have spectra with 
r ~ 2 or steeper at these energies. Conversely, where 
PGP emission is not strong, as in the GeV range, the 
latter processes may still prevail. For example, if the 
electron injection efficiency relative to the accretion ki- 
netic energy is f e = 0.01, the primary IC luminosity 
Lie - (21n7 e , max )" 1 / eJ L acc ~ 8 x 10 42 erg/s for our 
fiducial cluster, where 7 eimax ~ 3 x 10 is the maximum 
electron Lorentz factor determined from the balance be- 
tween acceleration and cooling times. This leads to flux 
~7x 10~ 12 erg cm -2 s~ 4 from D = 100 Mpc, close to 
the ~ 100 MeV EGRET upper limit for Coma (Reimer et 
al. 2003); nevertheless, PGP emission can be dominant 
at TeV if B s < 0.3£iG. Note that although f e is quite un- 
certain and values up to ~ 0.01 may be allowed, it may 
in fact be significantly lower. The spatial distribution of 
PGP emission should be similar to primary IC, track- 
ing the projected ring-like morphology of the accretion 
shock, so that p-p ir° gamma-rays may still be prominent 
for the central core regions in spatially resolved maps. 

To assess the detectability of PGP emission with var- 
ious observational facilities, we note that even though 
the full angular extent of the accretion shock diameter 
can be as large as 3.6°(i? s /3.2Mpc)(Z?/100Mpc)" 1 , the 
expected ring-like distribution is highly nonuniform; we 
thus choose to compare model fluxes with instrumental 
sensitivities for sources uniformly extended by 1°. Most 
promising for detection is the PGP IC component in TeV 
gamma-rays, where the currently operating HESS 4 array 
of Cerenkov telescopes can achieve point source sensitiv- 
ities ~ 10 _13 erg cm' 2 s _1 at 1 TeV in ^100 hr exposure 
time. With its 0.1° angular resolution, the 1° extended 
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source sensitivity should be roughly 10 _12 erg cm -2 s _1 , 
so that even our least luminous case of B s = ljuG may 
be detectable, not to mention the more luminous, lower 
B s cases. Furthermore, its 5° field of view may allow 
clear imaging of the annular morphology with a sin- 
gle pointing. Detecting the weaker PGP IC emission 
at GeV energies may be difficult even for the GLAST 5 
mission, with expected point source sensitivities reach- 
ing ~ 3 x 10~ 13 erg cm -2 s _1 at 1 GeV during a 1 yr 
survey. Moreover, it could likely be that other emission 
components such as primary IC show up above the PGP 
emission levels in this band (e.g. Miniati 2003; Gabici & 
Blasi 2004a). 

The PGP synchrotron component should be interesting 
for future facilities such as NeXT (Takahashi 2005), ca- 
pable of imaging hard X-rays in the 8-80 keV range at ~ 
30" resolution. The projected sensitivity for a 1° source 
is ~ 10 _12 erg cm -2 s _1 at 30 keV in a 100 ks exposure, 
so the cases of B s > 0.3/xG should be clearly detectable, 
and the accretion shock morphology may be studied in 
fine detail. The non-imaging HXD instrument on the 
soon-to-be launched ASTRO-E2 (Takahashi 2005), with 
a sensitivity of < 3 x 10 _12 erg cm -2 s" 1 at 20 keV, as 
well as instruments onboard INTEGRAL may also have 
chances of detection if B s > 1/iG. 

Finally, it is interesting to consider whether the model 
can explain the reported excess UV and hard X-ray emis- 
sion for the Coma cluster. In Fig|2| multifrcquency data 
compiled in Reimer et al. (2004) are compared with a 
model with B s = 0.7/iG, L p = 9x 10 45 erg/s and T inj = 4 
Gyr, otherwise being the same as in Fig^ Despite ap- 
pearing rather consistent, we recall that the UV and 
hard X-ray data correspond only to the inner < 1° re- 
gion, whereas the model emission should be extended 
out to projected radius ~ 1.8°. To explain the obser- 
vations with just the fraction of emission projected onto 
the core, more extreme parameters for L p and/or td yn 
may be necessary. Moreover, Bowyer et al. (2004) argue 
that the UV emission is spatially correlated with ther- 
mal X-rays, favoring a p-p secondary origin. In any case, 
further observations at these as well as gamma-ray ener- 
gies should provide conclusive answers. Note that neither 
the PGP nor p-p secondary synchrotron component from 
the accretion shock can account for the radio emission 
(Thierbach, Klein & Wielebinski 2003), calling for a dif- 
ferent origin such as particles accelerated by turbulence 
(Brunetti 2004; Fujita, Takizawa & Sarazin 2003). 

4. DISCUSSION 

Successful observations of the PGP emission discussed 
above entails a number of important implications, by pro- 
viding: 1) a clear signature of UHE protons in cluster ac- 
cretion shocks, and hence a test of particle acceleration 
theory on the largest scales; 2) a sensitive probe of mag- 
netic fields in the outermost regions of clusters, which is 
crucial for understanding the origin of intergalactic and 
intracluster magnetic fields, propagation of UHE cosmic 
rays, etc.; 3) a potential probe of the accretion shock it- 
self, which still lacks firm observational evidence despite 
being a robustly predicted cosmological phenomenon; 4) 
a useful tool for gamma-ray absorption studies of the 
IRB, owing to the hard, steady TeV spectra, as opposed 
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to the highly variable spectra of blazars. 

Our assumption of r\ = 1 as inferred for SNRs may 
not be completely justified for cluster accretion shocks. 
Larger r\ would decrease £ max and hence the amount 
of pair production and PGP emission; e.g. for the case 
of B s = 0.3/xG, 1) — 3 and 10 lead to IC luminosities 
lower by factors of ~ 3 and ~ 20, respectively. On the 
other hand, the proton injection efficiency can be higher 
than our choice of f p = 0.1. Increased emission may also 
result from proton spectra harder than 2 due to non- 
linear acceleration effects (Gabici & Blasi 2004b; Kang 
& Jones 2005), as well as from tertiary pairs produced 
in 7-7 interactions of PGP gamma-rays with the CMB 
(Timokhin, Aharonian & Neronov 2004). 

The present study may be more appropriate for clus- 
ters that are relatively relaxed rathen than those associ- 



ated with a current or recent major merger, although 
accretion shocks should be present even in the latter 
case (Ryu et al. 2003). The accretion process in reality 
is neither spherical nor steady, involving clumpy merg- 
ing along filamentary structures. However, as long as 
the merging leads to strong, high Mach number shocks, 
proton acceleration may proceed even through succes- 
sive merging events (Kang & Jones 2005). At any rate, 
more detailed investigations of the processes discussed 
here should be carried out by incorporating them in nu- 
merical simulations of large scale structure formation. 



We acknowledge fruitful discussions with E. Churazov, 
B. Dingus, T. Enfilin, S. Gabici, J. Kirk, F. Miniati, K. 
Nakazawa and T. Takahashi. 
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Fig. 1. — Spectra of proton-induced emission from our fiducial cluster accretion shock, for B a = 0.1, 0.3 and l^tG. The p-p 7r° component 
has been multiplied by 10. 
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Fig. 2. — Spectra of proton-induced emission for the model parameters labelled in the figure (see also text), compared with observational 
data for the Coma cluster from EUVE, BeppoSax, RXTE and EGRET compiled in Reimer et al. (2004). 



